Updating object motion dynamics from multiple radio signals

ABSTRACT

In one embodiment, a service receives signal characteristic data indicative of characteristics of wireless signals received by one or more antennas located in a particular area. The service uses the received signal characteristic data as input to a Bayesian inference model to predict physical states of an object located in the particular area. A physical state of the object is indicative of at least one of: a mass, a velocity, an acceleration, a surface area, or a location of the object. The service updates the Bayesian inference model based in part on the predicted state of the object and a change in the received signal characteristic data and based in part by enforcing Newtonian motion dynamics on the predicted physical states.

RELATED APPLICATION

This application claims priority to U.S. Provisional Patent Application No. 62/857,940, filed on Jun. 6, 2019, entitled “UPDATING OBJECT MOTION DYNAMICS FROM MULTIPLE RADIO SIGNALS” by Maluf et al., the contents of which are incorporated by reference herein.

TECHNICAL FIELD

The present disclosure relates generally to computer networks, and, more particularly, to updating object motion dynamics from multiple radio signals.

BACKGROUND

The problem of inferring an object model and the behavior of an object with a high resolution is not an easy task when doing so must be accomplished in an efficient, real time, and comprehensive manner. Indeed, such an undertaking requires near optimal integration of all of the available data from any number of data sources. More specifically, with the proliferation of wireless communications throughout the world, it now becomes possible to leverage the radio frequency (RF) communication signals to infer a model of the corresponding object. For example, a slight variance in the phase shift of a wireless signal, with no change in the attenuation of the signal, can be used to infer that the object is in some sort of motion. However, challenges still remain to translate the wireless signals into a usable object model that can be used for purposes of object recognition, delivering services to the object, etc.

BRIEF DESCRIPTION OF THE DRAWINGS

The embodiments herein may be better understood by referring to the following description in conjunction with the accompanying drawings in which like reference numerals indicate identically or functionally similar elements, of which:

FIG. 1 illustrate an example computer network;

FIG. 2 illustrates an example network device/node;

FIG. 3 illustrates an example wireless network;

FIG. 4 illustrates an example of representing a human as a center of mass;

FIG. 5 illustrates an example of representing a human using kinematics;

FIG. 6 illustrates an example of transmitting wireless signals towards a human;

FIGS. 7A-7C illustrate examples of wireless signals;

FIG. 8 illustrates an example Bayesian inference framework;

FIGS. 9A-9C illustrate example plots of test results; and

FIG. 10 illustrates an example simplified procedure for updating object motion dynamics based on the characteristics of received wireless signals.

DESCRIPTION OF EXAMPLE EMBODIMENTS Overview

According to one or more embodiments of the disclosure, a service receives signal characteristic data indicative of characteristics of wireless signals received by one or more antennas located in a particular area. The service uses the received signal characteristic data as input to a Bayesian inference model to predict physical states of an object located in the particular area. A physical state of the object is indicative of at least one of: a mass, a velocity, an acceleration, a surface area, or a location of the object. The service updates the Bayesian inference model based in part on the predicted state of the object and a change in the received signal characteristic data and based in part by enforcing Newtonian motion dynamics on the predicted physical states.

DESCRIPTION

A computer network is a geographically distributed collection of nodes interconnected by communication links and segments for transporting data between end nodes, such as personal computers and workstations, or other devices, such as sensors, etc. Many types of networks are available, ranging from local area networks (LANs) to wide area networks (WANs). LANs typically connect the nodes over dedicated private communications links located in the same general physical location, such as a building or campus. WANs, on the other hand, typically connect geographically dispersed nodes over long-distance communications links, such as common carrier telephone lines, optical lightpaths, synchronous optical networks (SONET), synchronous digital hierarchy (SDH) links, or Powerline Communications (PLC), and others. Other types of networks, such as field area networks (FANs), neighborhood area networks (NANs), personal area networks (PANs), etc. may also make up the components of any given computer network.

In various embodiments, computer networks may include an Internet of Things network. Loosely, the term “Internet of Things” or “IoT” (or “Internet of Everything” or “IoE”) refers to uniquely identifiable objects (things) and their virtual representations in a network-based architecture. In particular, the IoT involves the ability to connect more than just computers and communications devices, but rather the ability to connect “objects” in general, such as lights, appliances, vehicles, heating, ventilating, and air-conditioning (HVAC), windows and window shades and blinds, doors, locks, etc. The “Internet of Things” thus generally refers to the interconnection of objects (e.g., smart objects), such as sensors and actuators, over a computer network (e.g., via IP), which may be the public Internet or a private network.

Often, IoT networks operate within a shared-media mesh networks, such as wireless or PLC networks, etc., and are often on what is referred to as Low-Power and Lossy Networks (LLNs), which are a class of network in which both the routers and their interconnect are constrained. That is, LLN devices/routers typically operate with constraints, e.g., processing power, memory, and/or energy (battery), and their interconnects are characterized by, illustratively, high loss rates, low data rates, and/or instability. IoT networks are comprised of anything from a few dozen to thousands or even millions of devices, and support point-to-point traffic (between devices inside the network), point-to-multipoint traffic (from a central control point such as a root node to a subset of devices inside the network), and multipoint-to-point traffic (from devices inside the network towards a central control point).

Fog computing is a distributed approach of cloud implementation that acts as an intermediate layer from local networks (e.g., IoT networks) to the cloud (e.g., centralized and/or shared resources, as will be understood by those skilled in the art). That is, generally, fog computing entails using devices at the network edge to provide application services, including computation, networking, and storage, to the local nodes in the network, in contrast to cloud-based approaches that rely on remote data centers/cloud environments for the services. To this end, a fog node is a functional node that is deployed close to fog endpoints to provide computing, storage, and networking resources and services. Multiple fog nodes organized or configured together form a fog system, to implement a particular solution. Fog nodes and fog systems can have the same or complementary capabilities, in various implementations. That is, each individual fog node does not have to implement the entire spectrum of capabilities. Instead, the fog capabilities may be distributed across multiple fog nodes and systems, which may collaborate to help each other to provide the desired services. In other words, a fog system can include any number of virtualized services and/or data stores that are spread across the distributed fog nodes. This may include a master-slave configuration, publish-subscribe configuration, or peer-to-peer configuration.

Low power and Lossy Networks (LLNs), e.g., certain sensor networks, may be used in a myriad of applications such as for “Smart Grid” and “Smart Cities.” A number of challenges in LLNs have been presented, such as:

1) Links are generally lossy, such that a Packet Delivery Rate/Ratio (PDR) can dramatically vary due to various sources of interferences, e.g., considerably affecting the bit error rate (BER);

2) Links are generally low bandwidth, such that control plane traffic must generally be bounded and negligible compared to the low rate data traffic;

3) There are a number of use cases that require specifying a set of link and node metrics, some of them being dynamic, thus requiring specific smoothing functions to avoid routing instability, considerably draining bandwidth and energy;

4) Constraint-routing may be required by some applications, e.g., to establish routing paths that will avoid non-encrypted links, nodes running low on energy, etc.;

5) Scale of the networks may become very large, e.g., on the order of several thousands to millions of nodes; and

6) Nodes may be constrained with a low memory, a reduced processing capability, a low power supply (e.g., battery).

In other words, LLNs are a class of network in which both the routers and their interconnect are constrained: LLN routers typically operate with constraints, e.g., processing power, memory, and/or energy (battery), and their interconnects are characterized by, illustratively, high loss rates, low data rates, and/or instability. LLNs are comprised of anything from a few dozen and up to thousands or even millions of LLN routers, and support point-to-point traffic (between devices inside the LLN), point-to-multipoint traffic (from a central control point to a subset of devices inside the LLN) and multipoint-to-point traffic (from devices inside the LLN towards a central control point).

An example implementation of LLNs is an “Internet of Things” network. Loosely, the term “Internet of Things” or “IoT” may be used by those in the art to refer to uniquely identifiable objects (things) and their virtual representations in a network-based architecture. In particular, the next frontier in the evolution of the Internet is the ability to connect more than just computers and communications devices, but rather the ability to connect “objects” in general, such as lights, appliances, vehicles, HVAC (heating, ventilating, and air-conditioning), windows and window shades and blinds, doors, locks, etc. The “Internet of Things” thus generally refers to the interconnection of objects (e.g., smart objects), such as sensors and actuators, over a computer network (e.g., IP), which may be the Public Internet or a private network. Such devices have been used in the industry for decades, usually in the form of non-IP or proprietary protocols that are connected to IP networks by way of protocol translation gateways. With the emergence of a myriad of applications, such as the smart grid advanced metering infrastructure (AMI), smart cities, and building and industrial automation, and cars (e.g., that can interconnect millions of objects for sensing things like power quality, tire pressure, and temperature and that can actuate engines and lights), it has been of the utmost importance to extend the IP protocol suite for these networks.

FIG. 1 is a schematic block diagram of an example simplified computer network 100 illustratively comprising nodes/devices at various levels of the network, interconnected by various methods of communication. For instance, the links may be wired links or shared media (e.g., wireless links, PLC links, etc.) where certain nodes, such as, e.g., routers, sensors, computers, etc., may be in communication with other devices, e.g., based on connectivity, distance, signal strength, current operational status, location, etc.

Specifically, as shown in the example network 100, three illustrative layers are shown, namely the cloud 110, fog 120, and IoT device 130. Illustratively, the cloud 110 may comprise general connectivity via the Internet 112, and may contain one or more datacenters 114 with one or more centralized servers 116 or other devices, as will be appreciated by those skilled in the art. Within the fog layer 120, various fog nodes/devices 122 (e.g., with fog modules, described below) may execute various fog computing resources on network edge devices, as opposed to datacenter/cloud-based servers or on the endpoint nodes 132 themselves of the IoT layer 130. Data packets (e.g., traffic and/or messages sent between the devices/nodes) may be exchanged among the nodes/devices of the computer network 100 using predefined network communication protocols such as certain known wired protocols, wireless protocols, PLC protocols, or other shared-media protocols where appropriate. In this context, a protocol consists of a set of rules defining how the nodes interact with each other.

Those skilled in the art will understand that any number of nodes, devices, links, etc. may be used in the computer network, and that the view shown herein is for simplicity. Also, those skilled in the art will further understand that while the network is shown in a certain orientation, the network 100 is merely an example illustration that is not meant to limit the disclosure.

Data packets (e.g., traffic and/or messages) may be exchanged among the nodes/devices of the computer network 100 using predefined network communication protocols such as certain known wired protocols, wireless protocols (e.g., IEEE Std. 802.15.4, Wi-Fi, Bluetooth®, DECT-Ultra Low Energy, LoRa, etc.), PLC protocols, or other shared-media protocols where appropriate. In this context, a protocol consists of a set of rules defining how the nodes interact with each other.

FIG. 2 is a schematic block diagram of an example node/device 200 that may be used with one or more embodiments described herein, e.g., as any of the nodes or devices shown in FIG. 1 above or described in further detail below. The device 200 may comprise one or more network interfaces 210 (e.g., wired, wireless, PLC, etc.), at least one processor 220, and a memory 240 interconnected by a system bus 250, as well as a power supply 260 (e.g., battery, plug-in, etc.).

Network interface(s) 210 include the mechanical, electrical, and signaling circuitry for communicating data over links coupled to the network. The network interfaces 210 may be configured to transmit and/or receive data using a variety of different communication protocols, such as TCP/IP, UDP, etc. Note that the device 200 may have multiple different types of network connections 210, e.g., wireless and wired/physical connections, and that the view herein is merely for illustration. Also, while the network interface 210 is shown separately from power supply 260, for PLC the network interface 210 may communicate through the power supply 260, or may be an integral component of the power supply. In some specific configurations the PLC signal may be coupled to the power line feeding into the power supply.

The memory 240 comprises a plurality of storage locations that are addressable by the processor 220 and the network interfaces 210 for storing software programs and data structures associated with the embodiments described herein. The processor 220 may comprise hardware elements or hardware logic adapted to execute the software programs and manipulate the data structures 245. An operating system 242, portions of which are typically resident in memory 240 and executed by the processor, functionally organizes the device by, among other things, invoking operations in support of software processes and/or services executing on the device. These software processes/services may comprise an illustrative super-resolution process 248, as described herein. Note that while process 248 is shown in centralized memory 240 alternative embodiments provide for the process to be specifically operated within the network interface(s) 210.

It will be apparent to those skilled in the art that other processor and memory types, including various computer-readable media, may be used to store and execute program instructions pertaining to the techniques described herein. Also, while the description illustrates various processes, it is expressly contemplated that various processes may be embodied as modules configured to operate in accordance with the techniques herein (e.g., according to the functionality of a similar process). Further, while the processes have been shown separately, those skilled in the art will appreciate that processes may be routines or modules within other processes.

FIG. 3 illustrates an example wireless network 300, according to various embodiments. Wireless network 300 may be deployed to a physical location, such as floor 302 shown, and may include various infrastructure devices. These infrastructure devices may include, for example, one or more access points (APs) 304 that provide wireless connectivity to the various wireless clients 306 distributed throughout the location. For illustrative purposes, APs 304 a-304 d and clients 306 a-306 i are depicted in FIG. 3. However, as would be appreciated, a wireless network deployment may include any number of APs and clients.

A network backbone 310 may interconnect APs 304 and provide a connection between APs 304 and any number of supervisory devices or services that provide control over APs 304. For example, as shown, a wireless LAN controller (WLC) 312 may control some or all of APs 304 a-304 d, by setting their control parameters (e.g., max number of attached clients, channels used, wireless modes, etc.). Another supervisory service that oversees wireless network 300 may be a monitoring and analytics service 314 that measures and monitors the performance of wireless network 300 and, if so configured, may also adjust the operation of wireless network 300 based on the monitored performance (e.g., via WLC 312, etc.). In further embodiments, as detailed below, monitoring and analytics service 314 may also be configured to perform object modeling from a motion and/or kinematics standpoint for purposes such as object detection/identification, object tracking, object behavioral predictions, alerting, and the like. Note that service 314 may be implemented directly on WLC 312 or may operate in conjunction therewith, in various implementations.

Network backbone 310 may further provide connectivity between the infrastructure of the local network and a larger network, such as the Internet, a Multiprotocol Label Switching (MPLS) network, or the like. Accordingly, WLC 312 and/or monitoring and analytics service 314 may be located on the same local network as APs 304 or, alternatively, may be located remotely, such as in a remote datacenter, in the cloud, etc. To provide such connectivity, network backbone 310 may include any number of wired connections (e.g., Ethernet, optical, etc.) and/or wireless connections (e.g., cellular, etc.), as well as any number of networking devices (e.g., routers, switches, etc.).

The types and configurations of clients 306 in network 300 can vary greatly, ranging from powerful computing devices to any number of different types of IoT nodes/devices. For example, clients 306 a-306 i may include, but are not limited to, wireless sensors, actuators, thermostats, relays, wearable electronics, and the like.

—Motion Dynamics and Inverse-Kinematics Super-Resolution with Radio Sensing—

As noted above, inferring an object model and the behavior of an object from wireless signals now becomes possible thanks to the proliferation of wireless communications. For example, in a simple case, radio signals received by an AP 304 in wireless network 300 can be used to detect when a client 302 is approaching it (e.g., for purposes of making a physical security assessment, etc.). According to various embodiments, Bayesian inference methods are introduced herein that allow a service (e.g., service 314) to infer the motion dynamics and/or inverse-kinematics of an object with super-resolution, based on received wireless signals.

To illustrate the super-resolution modeling techniques introduced herein, consider the operation of the Global Positioning System (GPS). As would be appreciated, the data captured by a GPS receiver comprises time echoes from in-sight satellites. However, these time records are useless on their own. When combined with longitude, latitude, and elevation as state estimates, though, GPS functions as one of the most precise location and navigation systems ever. To perform the estimations in GPS, a Doppler Differential Kalman Filter is used, which estimates the probability distribution of interest for the current state conditioned on the readings/measurements taken prior to the current time.

The techniques herein, on the other hand, are directed to providing an optimal statistical inference defined by the physics of motion as the underlying model/prior merging the data produced by multiple streams (e.g., wireless signals from any number of devices), to derive new and unmeasured states of the modeled object characteristics. In particular, the present disclosure describes a number of techniques that apply a Bayesian inference method of super-resolution to infer: a.) motion dynamics of objects and/or b.) inverse kinematics of objects, where the inferences are based on the difference of the models' state estimations and the measured signals of any combination of either single-spectral radio signal data or multi-spectral radio signal data. Notably, the produced inference is deemed optimal, and therefore super-resolved for the motion. Such techniques are widely applicable, and may be particularly well suited for health care, security, and safety, among many other use-cases.

Specifically, in the context of modern radio imaging, the proliferation of radio technology and devices in consumers, industry and governments present a larger number and great diversity in radio emitters and sensors. The techniques herein may be applied to this growing paradigm of radio emitter and sensors through the dynamics and inverse kinematics state estimations of the physical states of the objects in question. For instance, in terms of motion dynamics, objects may be modeled to their intrinsic elements using the Newtonian physics of motion (motion dynamics). By way of example, FIG. 4 illustrates an example human 400 that can be represented by his or her center of mass (CM) 402 (e.g., center of gravity) on which Newtonian forces apply.

The kinematics behavior of the object can similarly be modeled by applying kinematics laws to the intrinsic element of the articulations (joints) of the object. The object's components are therefore assumed to have a rigid shape. For example, FIG. 5 illustrates a kinematics state representation of a person 500. As shown, person 500 may be represented as a series of chained points connected by rigid segments such as head 502, right (R.) shoulder 504, right elbow 506, right wrist 508, left (L.) shoulder 510, left elbow 512, left wrist 514, chest 516, abdomen 518, right hip 520, right knee 522, right ankle 524, left hip 526, left knee 528, and left ankle 530. As would be appreciated, the number and specificity of these points can vary, as desired. Each such point 502-530 may have its own six degrees of freedom in the kinematics model, but each segment may be constrained by the others. For example, the movement of right knee 522 may, in turn, cause right ankle 524 to also move. Thus, the movement of points 502-530 may be interdependent on one another.

Therefore, an example use case of the techniques herein would be to infer a human's motion dynamics (location) and inverse kinematics (physical behavior) from received radio signals. This particular example, for instance, specifies the human constitute of multiple points joints (kinematics) following the general motion dynamics. To illustrate the relationship between the motion dynamics and kinematics behavior of an object, assume that the service uses radio signals to determine that right hip 520 is moving in a specific direction. In turn, the service can leverage inverse kinematics to propagate this motion to the other points 502-518 and 522-530 in a constrained and statistical sense.

Said differently, the techniques herein allow for a service, such as service 312, to represent the motion (motion dynamics) and behavior (kinematics) operating upon data derived from radio imaging using super-resolution Bayesian methods and signal processing methods to model the electromagnetic field variability as it interacts with the objects. Though one possible outcome of these techniques is an avatar mimicking a remote person in presence of radio signals, other outcomes may be made to model objects using motion dynamics and inverse kinematics, accordingly. Note also that Digital Signal Processing referred to herein pertains to the nature of the raw signal measured at the antenna capable of detecting the signal amplitudes and the phase shifts of the electromagnetic signals.

By way of example, FIG. 6 illustrates a simplified example of APs 304 j and 304 k transmitting wireless signals 602 j and 602 k, respectively, towards person 500. In turn, each AP 304 may receive wireless signals in response, such as reflected signals, which the service associated with APs 304 can use for purposes of motion and inverse kinematics analysis.

The problem of integrating information from different radio sensors to answer particular questions is a familiar one in remote sensing, and often is referred to as “data fusion.” “Data fusion,” as commonly defined in the state of the art, however, is fundamentally incorrect as data should never be tampered with, let alone “fused.” Data is what was observed, and as such cannot be changed after the event. In contrast, the techniques herein operate under the theory that it is possible to construct physical, mathematical, and behavioral models (motion and kinematics models) of the underlying physics of an object that “best” predicts the observed data.

From such physical models, it is possible to project new derived data points at a higher resolution and fidelity, which appears as a set of systems' state variables. This set is a union of state variables that are measured and that are purely derived from physical and behavioral models. That is, given the model, one can probabilistically predict what would be observed (the expected data) and the non-observed through the underlying models.

In general, super-resolution is an optimal method that combines the sensed data with the physics models' optimal precisions. Continuing the GPS example above, by measuring the time between the GPS receiver and the satellites, GPS produces derived states, namely longitude, latitude, and elevation. Therefore, advanced mathematical methods are used in GPS to combine the data from multiple satellites as well the data accumulated over time to achieve high precision. In this case, combining the computed GPS state with the behavior of the physical object in terms of its motion as a function of time resolves further the actual position as the motion characteristics (e.g., acceleration, velocity, and displacement) and in this case the limits of laws of physics that bind the underlying moving object.

The extension of this analogy to motion dynamics and inverse kinematics would be a person equipped with multiple GPS receivers located at different points of their body, which would resolve further through kinematics models of the body and better precisions of the locations of the joints. Inverse kinematics super-resolution is a complex computation process depending on the motion and kinematics models acting as a strong prior. (Notably, body models are often defined statically by the joint components.)

The implication of a strong model/prior optimally eliminates noise from data. This is also known as the Kalman Filter. Therefore, for the analogy example according to the techniques herein, if the person is equipped two low resolution GPS receivers at different joint location, the inference links the two GPS data through inverse kinematics and would produce a higher resolution when compared to each independent GPS findings.

According to various embodiments, the proposed techniques herein utilize radio signals to infer positions of the articulation points of a modeled kinematics chained object. According to the techniques herein, an object may be super-resolved independently in time on its own motion as a function. The object motion is therefore inferred and super-resolved over time as a single dense object represented by its center mass location the and the object's mass. In addition, according to the techniques herein, the object kinematics may be inferred and further super-resolved over time based on the inverse-kinematics as defined by the chained objects. Further inference on the expected data computed from kinematics can be used to update further the model. For example, human motion kinematics has a known behavioral pattern (e.g., walking, sitting, etc.). The same inference, and when it differs from known kinematic or established behaviors, may constitute a clue of the behavior, to determine whether the model is trustworthy. For example, sudden acceleration of the object in the vertical direction may be presumed to be a fall, as it does not follow the known pattern.

The motion dynamics and kinematics models act as the central repository of all the radio information in the data and is constructed from the prior knowledge of physical processes and how physical dynamics interacts with their environment. This central model could loosely be described as the system that “fuses” data, but it is not itself “data” or a “data product”. The Bayesian probabilistic estimation approach used herein not only allows estimation of the most probable model given prior domain knowledge and the data, but it also estimates the uncertainty associated with the model. In particular, if the model's uncertainty is high, it means that there is insufficient data/prior knowledge to accurately construct the required model.

The Bayesian model-based data integration outlined above (a.k.a. Dynamic and Extended Kalman Filter), in principle solves the problem of how to integrate information from multiple sensors, namely radio antennas including multi-spectral and other sensing modality.

In multi-spectral sensor data, different resolutions in wavelength will update the models at different resolutions, enhancing further the model resolution. In this case, enhancing the resolution beyond the Nyquist limit of sensor resolution is what is referred to herein as “super-resolution.” In other words, super-resolution combines optimally the input data where the derived resolution exceeds the resolution of the highest sensor resolution.

Loosely said, super-resolution for motion dynamics and inverse-kinematics breaks the Nyquist limit. Single sensor super-resolution methods would hope to improve the resolution of the underlying data and overcoming the Nyquist limit only when motion is established. However, when more sensors and data points translated through their underlying physics and behavioral laws to the high-resolution models, high-resolution data is consequently produced with optimal statistical definition. Simply put, super-resolution is a tradeoff between the highly correlated statistical data against multi-dimensional interpolation and extrapolation through the underlying physical and behavioral models.

It is also quite useful for the service to identify changes when the difference between what the model behavior should be (estimating motion and kinematics) and the observed data is exceeding the noise threshold. The models can also be used in a differential approach against the data with the options to either resolve further additional models and new models' parameters or resolve and infer a larger parameter set (e.g., increase the number of joints for a complex object).

According to the illustrative embodiments herein, certain mathematics and formulas have been conceived for techniques herein, which are illustratively described below in the following Parts. In practice, uncertainty will be dominant in electromagnetic sensing, the choice of models and derivation are intended to carry forward the uncertainty and approximation errors over to the target medium models' characteristics. These uncertainties will skew the medium inferred expected parameters. Since the medium is a physical object, its parametric representation will be of Gaussian with a mean and variance. However, the skewedness of the uncertainty is readily removed from the final statistics as the model's nature are known a-priori. For example, humans may be the actual object medium in question and not a random object and therefore the model/prior is known to a certainty of fidelity for further inference in the inference derivation.

Notably, before continuing further in the description, it is important that the following symbols and nomenclature are understood:

-   -   1. Angular frequency ω (radians/sec)=2πf, where f=frequency (Hz)     -   2. This fall-off is far more rapid than the classical radiated         far-field E and B fields, which are proportional to the         inverse-distance (1/d) where d is the Euclidean distance.     -   3. The power density Γ (Gamma) of an electromagnetic wave is         proportional to the inverse of the square of distance from         source (Γ∝1/d²).

Part I: Signal Reflection and Surface Models

When an electromagnetic field is incident upon a boundary, in general it will split up into reflected and refracted fields. For a wave striking a separation interface of materials, Snell's law state that the incident angle equals the reflection angle θ1=02 such that: ν₁ sin θ₁=ν₂ sin θ₂ with ν₁, ν₂ being the refractive indices of the two media, the reflection coefficient is:

$R_{E} = {\frac{{v_{1}\cos\;\theta_{1}} - \sqrt{v_{2}^{2} - {v_{1}^{2}\sin^{2}\theta_{1}}}}{{v_{1}\cos\;\theta_{1}} + \sqrt{v_{2}^{2} - {v_{1}^{2}\sin^{2}\theta_{1}}}}.}$

It is of interest to model the specular reflectance when the electromagnetic field reflects over a convex closed boundary object. When the total surface of the object represents curves represented in spherical terms, the closed integral for R_(E) over the convex surface converges to a Lambertian Reflectance approximation.

Lambertian reflectance is the property that defines an ideal “matte” or diffusely reflecting surface. The apparent brightness of a Lambertian surface to a receiver antenna is the same regardless of the antennas angle. More technically, the surface's reflectance is isotropic, and the intensity obeys Lambert's cosine law.

The Lambertian approximation is a model of the perceived reflectance of a surface or a component of the wave intensity in a way that is as close to reality. The Lambertian model is further validated where the closed integral of the reflections in indoors environment includes the line-of-sight and multipaths reflections.

With further considerations that the natural random aspects and variations of the underlying surfaces. The assumption that the electromagnetic radio reflectance can be approximated with the Lambertian's reflectance surfaces is valid choice. For example, the assumption is true for general human surfaces.

The reflection of the surface of a medium over convex surface is: R _(E)=√{square root over (A _(e))}Φ where Φ_(m) is infinitesimal reflected radiation flux for the Lambertian model for a point surface such that:

$\Phi = {\sqrt{\rho}\frac{1}{R_{s}}\left( {{E_{s}\cos\alpha^{s}} + E_{A}} \right)\cos\alpha^{r}}$ where ρ is the albedo (material reflection coefficient between 0-1) of the surface, E_(s) is the perceived radiation intensity direct from the source antenna, E_(A) is the perceived ambient radiation (i.e. multipaths), and R is the distance driven diffusion factor from source antenna. E_(A) is also inversely proportional dependent on distance from the source (for non-infinite distance sources). Further, α^(s) and α^(r) are the angle from the surface norm towards the source and receiver respectively.

A general derivation would consist of the source and receiver having different locations. For this example, it can be assumed that both source and receiver are collocated. The close integral for the line-of-sight and multipaths flux is approximated as follows:

$\Phi = {\sqrt{\rho}\frac{1}{R}\left( {E_{s} + E_{A}} \right)}$

To compute the Lambertian flux bidirectional reflectance distribution function (BRDF) for an object, the projected A_(e) is defined as the effective surface area, which may be modeled experimentally, being half (½) of the total surface for the object. For volumetric computations, an assumption can be made that objects present complex curvatures which can be further approximated in spherical and/or cylindrical intrinsic forms. For example, for a fractional cylindrical form as an assumption, the approximation with a radius r and height and h is: A _(e) =rh∫ ₀ ^(π)sin(x)dx=2rh where r and h are considered random variables and have a Gaussians in nature. When motions occur, objects will produce change in their displacements of multiple random surfaces which are perceived by a receiving antenna. The receiving antenna can therefore sense an electromagnetic signal intensity change E_(R) being the sum of the contributions of moving surfaces represented by their respective effective areas and flux. The sum of the projections is therefore:

A_(e) = rh∫₀^(π)sin (x)dx = 2  rh and $R_{E} = {\sum\limits_{i}^{N}{\sqrt{A_{ei}}\Phi_{i}}}$

With the assumption that not everything moves instantaneously, the measurement of the signals at very high sampling rate (e.g. micro-seconds range) will isolate each motion independently as a micro-motion and for a virtual micro-surface in the Euclidean space with a respective area A_(e) and micro-displacement dλ at that moment, the displaced volume is therefore: V=A _(e) dλ

The micro-mass can be also computed for an assumed mass density. Other curvatures and shapes models can be considered such as spheres, ellipsoids instead of a cylindrical form. The mathematical derivation will follow a similar derivation to a cylinder.

By way of example, consider the case illustrated in FIG. 7A. As shown, assume that the object has a three dimensional (3D) micro-surface 702 that has an area greater than λ. During operation, both line-of-sight waves 704 and multipath waves 706 will reflect off of micro-surface 702 and may be received by an AP, such as AP 304 j shown. From this, the motion 708 (e.g., the micro-displacement dλ) of micro-surface 702 can be computed.

Part II: Signal Attenuations

In terms of the parameters of an electromagnetic wave and the medium through which it travels, the wave impedance may first be computed. The surface (entry) to surface (exit) wave impedance Z of the medium is defined as the ratio of electrical and magnetic field and is evaluated as:

$Z = {\frac{E_{A}}{H_{A}} = \sqrt{\frac{j\;\omega\;\mu}{\sigma + {j\;\omega\; ɛ}}}}$ where E_(A) is the electric field and H_(A) is the magnetic field for the ambient radiation and have complex representation. The impedance also is a complex number where:

-   -   μ is the magnetic permeability,     -   ε is the (real) electric permittivity,     -   σ is the electrical conductivity of the material the wave is         travelling through (corresponding to the imaginary component of         the permittivity multiplied by omega),     -   j is the imaginary unit, and     -   ω is the angular frequency of the wave.

The electromagnetic field in the far-field is independent of the antennas (e.g., APs 304, etc.). Since wave impedance is the ratio of the strength of the electric and magnetic fields, which in the far field are in phase with each other. Thus, the far-field impedance of free space is purely resistive and is given as:

$Z_{0} = {\frac{E}{B} = {\sqrt{\frac{\mu_{0}}{ɛ_{0}}} = \frac{1}{ɛ_{0}c}}}$

For the virtual medium positioned in the free space between the source and receiver, its impedance is a function of frequency. This case is more prominent in outdoor where the effect of multipath is minimal. The total impedance is therefore the total sum of impedances and:

$Z = {\frac{1}{ɛ_{0}c} + \sqrt{\frac{j\;\omega\;\mu}{\sigma + {j\;\omega\; ɛ}}}}$

When the medium is in motion, the receiver is assumed to measure conceptually a varying impedance, assuming the medium totally lies within the source and receiver antennas.

Part III: Measurements

Assume now that the transmitter is a uniform source and that the signal is reflected via the object micro-surface, follows Lambertian reflectance, and then arrives at the receiving antenna, as in the case in FIG. 7A. In such a case, the Lambertian reflectance model derived previously for indoors may be used to model this reflectance, although other derivations are also contemplated for other use cases. In terms of power representation in a free space, attenuation (e.g., the received power at the receiving antenna) is given by:

$\begin{matrix} {{P_{RX} = {\frac{A_{RX}}{4\pi R^{2}}P_{TX}}},} & (a) \end{matrix}$ where P_(TX) is the power reflected (transmitted) by the surface, and A_(RX) is the effective area of the receiving antenna and R is the distance between that object (surface) and the receiving end. In terms of electric field strength, this yields the following:

$\begin{matrix} {{P_{RX} = {\frac{{E_{r}}^{2}}{Z_{0}}A_{RX}}},{and}} & (b) \\ {{P_{TX} = {\frac{{E_{s}}^{2}}{Z_{0}}\rho\; A_{e}}},} & (c) \end{matrix}$ where |E_(r)| and |E_(s)| are the magnitudes of the electric field strength at the receiving and transmitting ends, respectively. Further, A_(e) denotes the effective area of the object, Z₀ is the impedance of free space, Z₀=120πΩ, and ρ is the albedo. In other words, ρA_(e) denotes the effective micro-surface area of the object. Substituting from Equations (b) and (c) into (a) yields:

${E_{R}} = {\frac{\sqrt{\rho\; A_{e}}}{2\sqrt{\pi}R}{{E_{S}}.}}$

Now, let M denote the normalized magnitude of electric field strengths as follows:

$M = {\frac{E_{R}}{E_{s}} = {\frac{\sqrt{\rho\; A_{e}}}{2\sqrt{\pi}R}.}}$

Furthermore, the phase difference between the received and transmitted electric field is defined as follows:

$\psi\overset{\Delta}{=}{\psi_{E_{r}} - {\psi_{E_{s}}.}}$

The measured phase is folded/wrapped due to the recurrence characteristic of phase. Therefore, in various embodiments, the service may transform the measured phase into the true, ‘unwrapped’ value. In one embodiment, the device may unwrap the measured phase using the following algorithm, which is represented in pseudocode:

Input: Measured (wrapped) phase ψ_(l) of L subcarriers 1 = 1: L. Output: Unwrapped phase ϕ₁ of L subcarriers; 1 = 1 : L. Initialize ϕ₁ = ψ₁ for 1 = 2 : L do if ψ_(l) − ψ_(l−1) > π then d = d + 1 end if ϕ_(l) − ψ_(l) − 2πd endfor

For example, in the case of Orthogonal Frequency Division Multiplexing (OFDM) subcarriers, the difference between the unwrapped phases can be represented as follows:

$\phi\overset{\Delta}{=}{{\phi_{E_{r}} - \phi_{E_{s}}} = {{{- \frac{2\pi}{\lambda}}2R} = {{- \frac{4\pi}{\lambda}}R}}}$ where λ is the wavelength such that

${\lambda = \frac{c}{f}},$ c is the speed of light, and f is the frequency of operation.

Part IV: Dynamic System Model and Motion

Once motion is achieved, changes in the signal strength/gain and its phase are perceived on the receiver antenna. When considering infinitesimal motion, the signal phase change (phase shift) is affected by the medium displacement travelled divided by the wavelength for normalization for the angular notation.

For example, as shown in plot 710 in FIG. 7B, assume that the receiving antenna receives both line of sight and multipath signals. As shown, signal C is leading signals D and E, signal D is leading signal E and lagging behind signal C, and signal E is lagging behind both signals C and D. In such a case, the service may determine both the magnitude and phase of each signal and add the signals, accordingly.

In FIG. 7C, when motion occurs, as represented by dλ, the service may determine the corresponding change in magnitude dM and change in phase dφ, in order to model the motion of the object. This can be performed, in some embodiments, on a per-antenna and/or per-device (AP) basis.

From the previous equations, the wave formation is a function of direct and multipath reflections with the effect of multipath traversing the medium is approximated as follows:

${E_{R}} = {\frac{\sqrt{\rho\; A_{e}}}{2\sqrt{\pi}R}{E_{s}}}$

Infinitesimally, the derivative matrix (Jacobian) across multiple receivers over the distances R and the area A_(e) where A′_(e)=ρA_(e), is therefore:

$\begin{matrix} {J = \left\lbrack {\frac{{dE}_{R}}{dR}\mspace{25mu}\frac{{dE}_{R}}{dA}} \right\rbrack} \\ {= {\begin{bmatrix} {- \frac{\sqrt{A_{e}^{\prime}}}{R_{1}^{2}}} & \frac{\sqrt{A_{e}^{\prime}}}{R_{1}R_{2}} & \ldots & \frac{\sqrt{A_{e}^{\prime}}}{R_{1}R_{N}} & \frac{1}{2R_{1}\sqrt{A_{e}^{\prime}}} \\ \frac{\sqrt{A_{e}^{\prime}}}{R_{2}R_{1}} & {- \frac{\sqrt{A_{e}^{\prime}}}{R_{2}^{2}}} & \ldots & \frac{\sqrt{A_{e}^{\prime}}}{R_{2}R_{N}} & \frac{1}{2R_{2}\sqrt{A_{e}^{\prime}}} \\ \vdots & \vdots & \ddots & \vdots & \vdots \\ \frac{\sqrt{A_{e}^{\prime}}}{R_{N}R_{1}} & \frac{\sqrt{A_{e}^{\prime}}}{R_{N}R_{2}} & \ldots & {- \frac{\sqrt{A_{e}^{\prime}}}{R_{N}^{2}}} & \frac{1}{2R_{N}\sqrt{A_{e}^{\prime}}} \end{bmatrix}.}} \end{matrix}$

One observation made herein is that the derivative is invariant of the object shape as long the object projected A_(e) sensitivity to motion is not marginal as to the effect of motion. Further derivation can be constructed to model A_(e) being affected by motion when the underlying medium is considered non-symmetrical in shape. Note that {right arrow over (μ)}, A_(e) and velocity {right arrow over (ν)} are the dynamic motion of a micro-object.

Part V: Minimization and Dynamic System Estimation from Multiple Antennas

The specific form of the general model parameters that are considered by the service may be parameterized by a set of 3-dimensional parameters vector {right arrow over (u)}(t) and medium micro-area A_(e)(t) and a velocity {right arrow over (v)}(t) in a three-dimensional space as follows: p({right arrow over (u)},{right arrow over (v)}A _(e) |E ₁(t) . . . E _(n)(t))∝p(E ₁ . . . E _(n) |{right arrow over (u)},{right arrow over (v)},A _(e))p({right arrow over (u)},{right arrow over (v)},A _(e)) where Ê₁({right arrow over (u)},{right arrow over (v)},A_(e)) denotes the synthesized signal from the model and σ_(e) ² is the noise variance. This gives way to a Bayesian inference model that can be solved using the following:

${p\left( {\left. {{E_{1}(t)}\mspace{14mu}\ldots\mspace{14mu}{E_{n}(t)}} \middle| \overset{\rightarrow}{u} \right.,\overset{\rightarrow}{y},A_{e}} \right)} \propto {\prod\limits_{i}{\exp\left( {- \frac{\sum\left( {E_{i} - {{\hat{E}}_{i}\left( {\overset{\rightarrow}{u},\overset{\rightarrow}{y},A_{e}} \right)}} \right)^{2}}{2\sigma_{e}^{2}}} \right.}}$ where Ê₁({right arrow over (u)},{right arrow over (v)},A_(e)) is estimated from the models stated above. In various embodiments, applying Bayes' theorem to the problem yields: p(R ₁ . . . R _(n) ,A _(e) |E ₁(t) . . . E _(n)(t))∝p(E ₁(t) . . . E _(n)(t)|R ₁ . . . R _(n) ,A _(e))p(R ₁ . . . R _(n) ,A _(e)) where R₁ . . . R_(n) are the distances between object and the antennas, A_(e) is the effective area, and E₁(t) is the measured signal at the receiving i^(th) antenna, as well as the following: p(x,y,z,s|R ₁ . . . R _(n) ,A _(e))∝(p(u|R ₁ . . . R _(n) ,A _(e))p(x,y,z) where x,y,z are the position of the surface (x, y, z) and are estimated by solving the following set of quadratic equations for x, y and z and as velocities, v_(x) _(i) , v_(x) _(i) , v_(z) _(i) as follows:

${{{\left( {x - x_{i}} \right)^{2} + \left( {y - y_{i}} \right)^{2} + \left( {z - z_{i}} \right)^{2}} = {sR_{i}^{2}}}{\left( v_{x_{i}} \right)^{2} + \left( v_{y_{i}} \right)^{2} + \left( v_{z_{i}} \right)^{2}}} = \left( \frac{\frac{\lambda}{4\pi}d\phi_{i}}{dt} \right)^{2}$ where {right arrow over (u)}(t) and {right arrow over (v)}(t) are the vector system dynamics of the location and velocity of the micro-surface, dϕ_(i) is the change in the phase measured at the i^(th) antenna, and s is a scale to accommodate for the unknown constants. This states that the posterior distribution of the motion dynamics parameters and the effective area are proportional to the likelihood (e.g., the probability of observing the data given the parameters) multiplied by the prior distribution on the motion dynamics parameters and effective area.

For illustrative purposes only, the single-spectrum derivations above are for any numbers of antennas. However, those skilled in the art will appreciate that multiple spectral antennas derivations may be combined in a similar fashion as an extension to the single-spectrum derivations illustrated herein, and without departing from the intent and scope of the techniques embodied within the present disclosure. For example, signal data received by 2.4 GHz and 5 GHz antennas can be combined to make a single inference for a particular object.

Updating Object Motion Dynamics from Multiple Radio Signals

The present disclosure is directed to certain aspects of the techniques described above, particularly in relation to updating object motion dynamics from multiple radio signals. Indeed, using the above derivations and modeling approaches detailed above, a service for a wireless network may infer and update the various motion dynamics of a

Illustratively, the techniques described herein may be performed by hardware, software, and/or firmware, such as in accordance with the illustrative super-resolution process 248, which may include computer executable instructions executed by the processor 220 (or independent processor of interfaces 210) to perform functions relating to the techniques described herein (e.g., to provide a super-resolution service to a network).

Specifically, according to various embodiments, a service receives signal characteristic data indicative of characteristics of wireless signals received by one or more antennas located in a particular area. The service uses the received signal characteristic data as input to a Bayesian inference model to predict physical states of an object located in the particular area. A physical state of the object is indicative of at least one of: a mass, a velocity, an acceleration, a surface area, or a location of the object. The service updates the Bayesian inference model based in part on the predicted state of the object and a change in the received signal characteristic data and based in part by enforcing Newtonian motion dynamics on the predicted physical states.

Operationally, in various embodiments, the above techniques can be augmented by the super-resolution service/process enforcing the laws of motion dynamics on the motion of a detected body/object. Indeed, when using a Bayesian inference model to infer the physical state(s) of an object based on wireless signals received by any number of antennas, the inferred physical state(s) should still conform to the laws of Newtonian dynamics. Accordingly, in various embodiments, the Bayesian inference model may also take into account these dynamics laws, when making its predictions. For example, in one embodiment, the computed error for the model may also be based on whether the predicted object state is consistent with these laws.

As would be appreciated, the dynamics laws that the service may enforce are a function of the specific state parameter(s) predicted by the inference model, such as the location of the object, the mass of the object, the velocity of the object, and/or the acceleration of the object. Non-limiting examples of these constrains may include, but are not limited to, any or all of the following:

TABLE 1 Predicted Object State Parameter Constraint Surface Area The predicted surface area must be consistent with the estimated distances of the object to the signal receivers (e.g., APs) and motion of the object. Location The predicted location must be consistent with the previously-determined location of the object and/or the other predicted state parameter(s). For example, once the object appears within the beam of an antenna, the object cannot simply appear or disappear. In addition, the existence of the object in the beam must obey an entry and exit point from the beam. The predicted location of the object should also be consistent with its previous location and its velocity and/or acceleration. Mass The predicted mass of the object must be consistent with any other predicted state parameters. For example, the linear momentum of the object (e.g., Mass × Velocity) should be consistent with the force acted on the object (e.g., Mass × Acceleration), etc. Velocity The predicted velocity of the object must also be consistent with any other predicted state parameters. For example, the predicted velocity should be consistent with the predicted change in location of the object over time, the predicted acceleration of the object, the kinetic energy of the object (e.g., Mass × Velocity), the direction of motion of the object, etc. Acceleration The predicted acceleration of the object should also be consistent with any other predicted state parameters, as well as any expected values. For example, if the object is falling, the predicted acceleration may be constrained to not exceed that of gravitational acceleration (e.g., 9.8 m/s²).

In addition to enforcing Newtonian motion dynamics on the predicted object states, the service may also resolve the estimates of each individual body/object using the motion dynamics and measured data computed from the change in gain in the antenna(s), the phase shift rate of change, and/or the ratio of the antenna change(s) in gain over the antenna phase shift rate(s) of change.

In one embodiment, the service may employ a Kalman filter to continuously update the model. For example, such a filter may update the model parameters (e.g., mass, location/position, velocity, acceleration, etc.) in real time by continuously solving for the error between the motion dynamics prediction and the change in gain in the antennas, the phase shift rate of change, and ratio of the antenna change in gain over the antenna phase shift rate of change.

FIG. 8 illustrates an example Bayesian framework 800 for inferring object dynamics from multiple radio signals, according to various embodiments. For example, super-resolution process 248 may provide a service to the network by implementing Bayesian framework 800 (e.g., a non-linear Kalman filter or the like), to infer the characteristics of a given object, such as its motion dynamics.

As shown, framework 800 may include a Bayesian model 806 of the form described in Part V above. During operation, a simulation 808 may be performed using the prior of model 806, to form synthetic data 810. Synthetic data 810 is compared with the measured data 802 (e.g., the characteristics of the received signals), to form an inverse estimate 804. In turn, inverse estimate 804 is used to update model 806, allowing for repeated convergence over time.

In a further embodiment, the service could also identify the motion dynamics of the object/body using a digital twin of the assumed object in the real world. For example, if the service determines that a human is likely present within the area of interest, based on the received wireless signals, the service could create a digital twin of the person and ensure that any predictions are consistent with the expected movement of the digital twin.

Two example embodiments for which the techniques herein may be used, among any number of others, are the deriving, from radio signals: a.) avatar kinematics articulations and visualizations, and b.) object behavior classifications for object kinematics.

Regarding the first example, avatar kinematics articulations and visualizations derived from radio signals, avatars are digital models with multiple possibilities of graphical representations of people and objects and also rendering at large-scale mimicking the physical reality. Avatars are known in computing and especially in gaming. Conceptually, avatar models are interactive to input that is either controlled through user interaction interfaces such as keyboards or advanced haptics, or automated through computer simulations and other reasoning capabilities.

The example of the present disclosure stipulates that the interfaces of the avatars accept input from digital models in 2-Dimensions and 3-Dimensions, with physical movement in the dimension of time. These models are to obey the law of physics dynamics of motion and kinematics of the linked components of the models.

The techniques herein further use the methods of Bayesian super-resolution inference to assert, create, and update the fidelity and resolution of the avatars' models in real time of the people or objects. Avatar resolution is enhanced overtime as the continuous real-world measurement are processed through radio frequency signals when combined with the intrinsic simulations of already inferred models. These simulations constitute physics simulation of a free body motion dynamics and kinematics once the model type (e.g., human, dog, etc.).

For example, once a body is assumed to be moving horizontally, the law of kinetics energy and motion dynamics are used to extrapolate the next instant which includes locations, velocities, and accelerations which are rendered through the avatar. It is understood that the mentioned estimated physical quantities are in plural, to reflect the estimated quantities of all linked components of a body (e.g., hand, foot, head, etc.).

The avatars reflect the super-resolved estimated physical quantities based on the optimal calculations taking in consideration extrapolation calculations (projection) with the estimation of the models' errors and the radio frequency (RF) sensing data with the estimation of sensing errors.

Avatars models follow known physical models with connected chains of components with known range of behavior. For example, if the avatar is a human type, the kinematics of the human body kinematics is known and established as a model/prior by the components characteristics as well as the limitation that the articulation points offer to the movable component relative to what it is connected to. For example, the kinematics of a moving finger has a known range, whereas the finger is connected to a wrist with a known range whereas the wrist is connected to forearm with known range, etc.

Super-resolution with an avatar in context takes two orthogonal directions. For example, assuming the avatar is human, sensing data over time can articulate the motion of the avatar in a direction such that additional the number of the articulations points of the avatar are identified over time as well as the physical dimensions of the components constituting the avatar. Once avatars are established and simulated at large scale such as an industrial environment, the models are continuously inferred and updated with the radio signals at hand. Avatars states may also be further resolved through computer simulation constrained by the physical setting they are located in.

Originally, an avatar is assigned a mass and volume with a wide error margin. with time, the estimation of the mass and volume become more precise for the rendition and simulation as well interaction with the environment. For example, if the setting is an industrial manufacturing setting and the avatars represent the workers, the motion of the avatars is computed against the constraining environment such as walking horizontally, cannot traverse walls but must enter through doors, etc. These constraints are important to the inference to further super-resolve and remove the errors made in the inference due to poor or noisy data.

The avatar interaction with the environment constitute a feedback to the inference from the stage setting. Any failure in the inference engine using data while interacting the models reflects two possibilities: Either the data is poor to use, or the models are not correct to start with. These scenarios are flagged for correction. Either better antennas are place or the environment models are corrected if defects are found.

The advantage of avatars kinematics, and therefore visualization, offer advanced safety to many industrial and health care verticals. Industrial safety, for example, may be used to safeguard people from hazardous equipment. Health care verticals offer the benefit of observing patients in hospital settings. (Many other use cases may be presented, and the techniques herein are not limited to those examples specifically pointed out herein.) Many existing solutions in safety and surveillance are invasive in nature. Camera systems are known to invade privacy. Reconstructed avatars, on the other hand, offer a simulated behavior deriving from radio signals protecting privacy to the extent of the object or the individual's physical characteristics such as height, body mass index, and an approximate body composition.

Regarding the second example, object behavior classifications for object kinematics derived from radio signals, object behavior classifications follows the object of kinematics modeling and simulation. Object behaviors are mostly known and there are always unknown patterns. For example, if the body is human and the kinematics motion shows a sudden drop in center of gravity of the human (sudden acceleration and deceleration) into a sudden stop, such behavior maybe be labeled as a “fall”. The variety of the behavior classifications are wide to the motion kinematics of the body limited to the resolution made available through the inference. These behaviors could be organized in any arbitrary taxonomies such as “exercising” being a superset of “running” being a super set of the center of gravity displacement larger than 2.5 mph in combination of periodic oscillation of the center of gravity in the vertical direction over 3 inches. Object behavior classifications are therefore pattern matching of known behavior against the kinematics of the inferred models.

A prototype system was constructed using the techniques herein and used to make various inferences regarding an object in a wireless network. More specifically, data was collected from a plurality of APs operating on the 5 GHz frequency band. The three APs were arranged in a triangle with the first AP assumed to be the origin. The second AP was aligned with the first AP in the y-direction and located two meters away from the first AP in the x-direction. Finally, the third AP was located two meters away from the first AP in both the x-direction and y-direction. The object/client was located within the formed triangle and pinged the APs at an interval with a mean of two milliseconds and a variance of one millisecond. The received channel state information (CSI) was then filtered based on the MAC address of the client. During testing, multiple people were moving casually around the three APs.

FIGS. 9A-9C illustrate various examples of inferences made about the client during testing, according to various embodiments. Plot 900 in FIG. 9A shows the estimated distances R₁, R₂, and R₃ from the three APs, respectively, to the effective moving surface of the client. Based on these distances and the computations detailed above, the x-y positions of the client were then inferred, as shown in plot 910 in FIG. 9B. Plot 920 in FIG. 9C shows the effective surface area of the moving surface versus the samples. As would be appreciated, the effective surface area reflects the size of the moving surface.

FIG. 10 illustrates an example simplified procedure for updating object motion dynamics based on the characteristics of received wireless signals, in accordance with one or more embodiments described herein. For example, a non-generic, specifically configured device (e.g., device 200) may perform procedure 1000 by executing stored instructions (e.g., process 248), to provide a service to one or more networks. The procedure 1000 may start at step 1005, and continues to step 1010, where, as described in greater detail above, the service may receive signal characteristic data indicative of characteristics of wireless signals received by one or more antennas located in a particular area. In one embodiment, the one or more antennas are wireless access point antennas in a wireless network. In a further embodiment, the wireless signals may be received by a plurality of different antennas. In yet another embodiment, the characteristics of the wireless signal comprise one or more of: a change in gain of the wireless signals, a phase shift rate of change of the wireless signals, or a ratio of the change in gain over the phase shift rate of change of the wireless signals.

At step 1015, as detailed above, the service may use the received signal characteristic data as input to a Bayesian inference model to predict physical states of an object located in the particular area. In various embodiments, the physical state of the object may be indicative of at least one of: a mass, a velocity, an acceleration, or a location of the object. Notably, a Bayesian inference model can be used to predict the dynamics of an object present within the beam(s) of any number of antennas with super resolution.

At step 1020, the service may update the Bayesian inference model based in part on the predicted state of the object and a change in the received signal characteristic data and based in part by enforcing Newtonian motion dynamics on the predicted physical states. In various embodiments, the service may do so using a Kalman filter to solve for an error between the predicted physical states of the object and the characteristics of the wireless signals. In general, the service may enforce the motion dynamics on the predicted physical states by ensuring that any predicted physical states of an object (e.g., its location, mass, velocity, acceleration, etc.) are consistent with one another and/or any prior states of the object. For example, the service may ensure that a predicted location of the object is consistent with its previous location(s), such as by ensuring that the object does not suddenly disappear, the object has an associated entry and/or exit point within the beam, etc. In another example, if the service predicts both the mass and acceleration of the object, the service may also ensure that these predictions are consistent with one another according to their Newtonian relationship, e.g., force=mass×acceleration. Procedure 1000 then ends at step 1025.

It should be noted that while certain steps within procedure 1000 may be optional as described above, the steps shown in FIG. 10 are merely examples for illustration, and certain other steps may be included or excluded as desired. Further, while a particular order of the steps is shown, this ordering is merely illustrative, and any suitable arrangement of the steps may be utilized without departing from the scope of the embodiments herein.

The techniques described herein, therefore, provide for updating object motion dynamics from multiple radio signals. Notably, the technology herein, as it relates generally to state estimation for motion dynamics and inverse kinematics, addresses the use of single or multi-spectral radio signal data integration for businesses such as healthcare (e.g., end-to-end monitoring in healthcare, such as for in-patient and outpatient “man-down” monitoring) as well as physical security (e.g., tracking of moving objects without camera-based privacy concerns). In addition, the techniques herein allow for the determination of micro-motions of an object, localizations of the object (e.g., the x-y coordinates of the object in the area), and clustering of the localizations for a given time window (e.g., to assess the spread function of the coordinates). Such computations can be used, for example, to generate 3-D or 2-D renderings of the object, forming a heatmap of potential locations of the object, and the like.

While there have been shown and described illustrative embodiments that are directed to particular embodiments and features of motion dynamics and inverse-kinematics super-resolution with radio sensing, it is to be understood that various other adaptations and modifications may be made within the intent and scope of the embodiments herein. For example, while certain implementations are disclosed herein, such as avatars, human body sensing, access point sensors, and so on, the techniques herein are not limited as such and can be applied to any number of different types of arrangements and/or verticals. Further, while certain specific formulas have been shown and described, such formulas may be modified or presented in a different manner, and the techniques herein are not so limited.

The foregoing description has been directed to specific embodiments. It will be apparent, however, that other variations and modifications may be made to the described embodiments, with the attainment of some or all of their advantages. For instance, it is expressly contemplated that the components and/or elements described herein can be implemented as software being stored on a tangible (non-transitory) computer-readable medium (e.g., disks/CDs/RAM/EEPROM/etc.) having program instructions executing on a computer, hardware, firmware, or a combination thereof. Accordingly, this description is to be taken only by way of example and not to otherwise limit the scope of the embodiments herein. Therefore, it is the object of the appended claims to cover all such variations and modifications as come within the true intent and scope of the embodiments herein. 

What is claimed is:
 1. A method comprising: receiving signal characteristic data indicative of characteristics of wireless signals received by one or more antennas located in a particular area; using the received signal characteristic data as input to a Bayesian inference model to predict physical states of an object located in the particular area, wherein a physical state of the object is indicative of at least one of: a mass, a velocity, an acceleration, a surface area, or a location of the object; updating the Bayesian inference model based in part on the predicted state of the object and a change in the received signal characteristic data and based in part by enforcing Newtonian motion dynamics on the predicted physical states.
 2. The method as in claim 1, wherein the one or more antennas are wireless access point antennas in a wireless network.
 3. The method as in claim 1, wherein the characteristics of the wireless signal comprise one or more of: a change in gain of the wireless signals, a phase shift rate of change of the wireless signals, or a ratio of the change in gain over the phase shift rate of change of the wireless signals.
 4. The method as in claim 1, wherein updating the Bayesian inference model further comprises: using a Kalman filter to solve for an error between the predicted physical states of the object and the characteristics of the wireless signals.
 5. The method as in claim 1, wherein enforcing the Newtonian motion dynamics on the predicted physical states comprises: ensuring that a predicted location of the object is consistent with a previous location of the object.
 6. The method as in claim 1, wherein enforcing the Newtonian motion dynamics on the predicted physical states comprises: ensuring that a predicted mass and a predicted acceleration of the object are consistent with one another.
 7. The method as in claim 1, wherein the wireless signals are received by a plurality of different antennas.
 8. The method as in claim 1, further comprising: determining and based on the predicted physical states of the object, that the object is falling.
 9. An apparatus, comprising: one or more network interfaces to communicate with one or more networks; a processor coupled to the network interfaces and configured to execute one or more processes; and a memory configured to store a process executable by the processor, the process when executed configured to: receive signal characteristic data indicative of characteristics of wireless signals received by one or more antennas located in a particular area; use the received signal characteristic data as input to a Bayesian inference model to predict physical states of an object located in the particular area, wherein a physical state of the object is indicative of at least one of: a mass, a velocity, an acceleration, a surface area, or a location of the object; update the Bayesian inference model based in part on the predicted state of the object and a change in the received signal characteristic data and based in part by enforcing Newtonian motion dynamics on the predicted physical states.
 10. The apparatus as in claim 9, wherein the one or more antennas are wireless access point antennas in a wireless network.
 11. The apparatus as in claim 9, wherein the characteristics of the wireless signal comprise one or more of: a change in gain of the wireless signals, a phase shift rate of change of the wireless signals, or a ratio of the change in gain over the phase shift rate of change of the wireless signals.
 12. The apparatus as in claim 9, wherein the apparatus updates the Bayesian inference model further by: using a Kalman filter to solve for an error between the predicted physical states of the object and the characteristics of the wireless signals.
 13. The apparatus as in claim 9, wherein the apparatus enforces the Newtonian motion dynamics on the predicted physical states by: ensuring that a predicted location of the object is consistent with a previous location of the object.
 14. The apparatus as in claim 9, wherein the apparatus enforces the Newtonian motion dynamics on the predicted physical states by: ensuring that a predicted mass and a predicted acceleration of the object are consistent with one another.
 15. The apparatus as in claim 9, wherein the wireless signals are received by a plurality of different antennas.
 16. The apparatus as in claim 9, wherein the process when executed is further configured to: determine, based on the predicted physical states of the object, that the object is falling.
 17. A tangible, non-transitory, computer-readable medium storing program instructions that cause a device to execute a process comprising: receiving signal characteristic data indicative of characteristics of wireless signals received by one or more antennas located in a particular area; using the received signal characteristic data as input to a Bayesian inference model to predict physical states of an object located in the particular area, wherein a physical state of the object is indicative of at least one of: a mass, a velocity, an acceleration, a surface area, or a location of the object; updating the Bayesian inference model based in part on the predicted state of the object and a change in the received signal characteristic data and based in part by enforcing Newtonian motion dynamics on the predicted physical states.
 18. The computer-readable medium as in claim 17, wherein the one or more antennas are wireless access point antennas in a wireless network.
 19. The computer-readable medium as in claim 17, wherein the characteristics of the wireless signal comprise one or more of: a change in gain of the wireless signals, a phase shift rate of change of the wireless signals, or a ratio of the change in gain over the phase shift rate of change of the wireless signals.
 20. The computer-readable medium as in claim 17, wherein updating the Bayesian inference model further comprises: using a Kalman filter to solve for an error between the predicted physical states of the object and the characteristics of the wireless signals. 